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Introduction 


The continuing and recently accentuated search 
for aerodynamic ally efficient supersonic transport 
designs necessitates not only the fullest use of 
existing technology, but also requires the develop- 
ment of new analytical methods of evaluating poten- 
tially efficient configurations. Because of the 
large portion of total airplane drag associated 
vita the generation of lift, there is need for a 
thorough study of design concepts offering the pos- 
sibility of reduction of this drag component. 


The wing planform, which is of primary impor- 
tance in its effect on drag at lifting conditions, 
in past studies has been largely restricted to 
straight line leading and trailing edges. Such 
limitations have resulted mainly from a lack of 
analytical methods for estimating the aerodynamic 
characteristics of wings of arbitrary planform. 

This deficiency may now be overcome by application 
of modern digit al computers to the solution of 
linearized theory integral equations for wing plan- 
forms which may employ curved or cranked leading 
and trailing edges. Computer programs such as 
those discussed in this paper greatly expand the 
possibilities for the development of truly effi- 
cient supersonic airplane configurations. 

(O c 7 This paper illustrates the application of 
numerical methods to wing camber surfaces of arbi- 
trary planform, to obtain the surface shape 
required to support a specified pressure distribu- 


tion or, inversely, to obtain the pressure distri- 
bution on wings of specified shape. Two special 
cases are noted and demonstrated: in the direct 

solution, the surface shape corresponding to an 
optimum combination of loadings (for least drag at 
a specified lift) may be obtained; in the inverse 
solution, the pressure distribution on a flat wing 
of arbitrary planform may be obtained. 


Used in combination, the numerical methods 
allow the determination of linear theory drag-due- 
to-lift polars and lift-moment relationships for 
wings of arbitrary planform which may employ any 
specified surface shape. With attention given to 
the real flow considerations that limit the appli- 
cation of linear theory, 1 studies of the aerody- 
namic characteristics of arbitrary wing planforms 
with various surface shapes may be analytically 
conducted. 

To illustrate the use of the method, a set of 
examples is presented for a typical planform 
series. The agreement obtained between the numer- 
ical method and experiment is shown in a second set 
of examples, employing both flat and 
surfaces. 
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Symbols 


a(l,n) 

b 

C D 

C D,min 

Cm 

Cm^F 

ACm = Cjjj 

°L 

c L,a 

ACp 

l 

L,N 

L*,N* 

M 

R 

x,y,z 

x cp 

z c 

Azp = z P 


leading-edge grid element weighting 
factor 

wing span 

drag coefficient 

zero -lift drag coefficient of flat -wing 
configuration 

drag coefficient due to lift, 

C D " c D,min 

pitching -moment coefficient 

zero-lift pitching -moment coefficient 
of flat -wing configuration 

Cm^F 

lift coefficient 

lift- curve slope per degree angle of 
attack 

pressure coefficient 

lifting pressure coefficient 

overall length of wing measured in 
streamwise direction 

designation of influencing grid 
elements 

designation of field-point grid 
elements 

Mach number 

average value of influence function 
within a grid element 

Cartesian coordinate system, x-axis 
streamwise 

x- coordinate of wing center of pressure 
camber surface ordinate 
z c,Ze 

wing angle of attack, deg 


p = /m 2 - 1 


|_,T] dummy variables of integration for x 

and y, respectively 

t designates a region of integration 

bounded by the wing planform and the 
fore Mach cone from the point x,y 

A wing leading-edge sweepback angle 

Subscripts: 


designates various drag' components 
(see fig. 7) 


Discussion 

A typical wing planform described by a rec- 
tangular Cartesian coordinate system is illustrated 
in figure 1. Overlaid upon the wing planform is 
the grid system used in the numerical solution of 
the linear theory integral equation. A mosaic of 
whole and partial, grid elements, approximates the 
actual wing planform and surface shape. In 
accordance with the assumptions of linear theory, 
the wing has negligible thickness, and lies essen- 
tially in the z = 0 plane. 

In the numerical approach, the grid elements, 
identified by L and N, are arranged such that 
L is numerically equal to x and N is numeri- 
cally equal to Py, where x and py take on only 
integer values. Partial grid elements along the 
wing leading and trailing edges are used to permit 
a closer approximation to the actual wing planform. 
The grid system of figure 1 is rather coarse for 
illustrative purposes; in actual usage, many more 
grid elements are employed. 

With respect to a specified field point x,y, 
the upstream region of influence t (bounded by 
the fore Mach cone from x,y and the wing leading 
edge) , is approximated by the shaded grid elements 
of figure 1. Each of these elements has associ- 
ated with it an influence factor R which relates 
the effect of the element and its average pressure 
to the surface slope required to obtain_a specified 
lifting pressure at x,y. The factor R is deter- 
mined from an approximate solution to the linear 
theory integrals over the region bounded by the 
individual elements. 

The variation of R within the fore Mach 
lines from the field point and the corresponding 
grid element (L*,N*) is illustrated in figure 2. 

For a set of elements at a constant L* - L value, 
the sum of the R values is zero, the single neg- 
ative value at N* = N balancing all the others. 
Where L* = L, only a single_element is contained 
within the Mach lines, and R = 0. 

The basic equation relating the required sur- 
face slope at x,y to a prescribed pressure dis- 
tribution, written in the form of the numerical 
solution, is equation (l) of figure 2. 2 The effect 
of the upstream elements within t is included in 


“ 

WF 

FF 


the summation term, which adds directly to the pre- 
scribed lifting pressure coefficient at the field 
element to define the necessary, surface slope. In 
this equation, the factor A(L,N) accounts for 
partial grid elements, being equal to the element 
length in the x-direction. 

Equation (2) , a rearrangement of equation (l) , 
provides a solution to the inverse problem, that of 
solving for the lifting pressure coefficient corre- 
sponding to a specified surface shape. 5 in this 
equation the summation term utilizes previously 
determined lifting pressure coefficients, which are 
computed following a prescribed order of calcula- 
tion, i.e., from apex aft. With this procedure, 
and utilizing the fact that R(N*,L*) is zero, no 
unknown pressure coefficients arise in the numeri- 
cal summations. 

The wing surface slope variation required in 
equation (2) is provided by supplemental calcula- 
tions, which consist of determining the streamwise 
inclination of all grid elements from a set of cam- 
ber surface ordinates. The flat wing at a small 
angle of attack represents a special case of equa- 

tion (2), In which is a constant, 

dx 

The precision of the numerical method in 
defining the theoretical surface shape required to 
support a specified pressure distribution on wings 
of arbitrary planform has been illustrated in an 
NASA report. 2 Similar illustrations involving the 
inverse solution are contained in an NASA prospec- 
tive report, 5 from which a typical example is 
included here to show the type of detailed pressure 
distributions obtained from the numerical method, 
and the departures from more rigorous solutions. 

In figure 3 numerical solution pressure distribu- 
tions for a flat double-delta planform at two Mach 
numbers are compared with results from a super- 
position analysis. ^ Agreement between the two 
pressure distributions is generally quite good, 
although the numerical solution does not produce 
the sharp pressure peaks along Mach lines that 
characterize the superposition analysis. Compari- 
sons of total wing lift coefficient and center of 
pressure between the two solutions show reasonably 
good agreement as follows: 



M = lAl4 



M = 1.667 


c L,tt 

x cp / * 

c L,a 

x cp /* 

Superposition 

analysis 

0 . 051 ^ 

0.682 
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Not 

given 

Numerical 
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The use of the numerical methods to obtain 
camber surfaces for wings having a prescribed 
loading, and pressure distributions for flat wings 
of the same planform is illustrated in figures 4 
and 5 for a series of wings having the same span 
and length. The three wings have delta, ogee, and 
blunt ogive planforms, oriented with respect to the 
apex Mach lines such that the sweepback parameter 
(p cot A) of the delta wing is 0.70. 

In figure 4, the surface shape and loading 
distributions are shown for an optimum combination 
(having least drag at a specified lift coefficient) 




of a uniform, linear spanwise, and linear chordwise 
loading, calculated by the numerical method. The 
corresponding drag-due-to-lift factors are also 
tabulated. The general characteristics of the 
required surface shape for the delta and ogee plan- 
forms are quite similar, both exhibiting the 
extreme root incidence predicted by linear theory 
for wings with a sharp apex. The ogee wing has 
more abrupt surface curvature near the wing root, 
however, and remains nearly flat over the outboard 
60-percent span. In contrast, the blunt ogive 
planfprm exhibits much more moderate surface slopes, 
and less variation in surface pressures. The drag- 
due-to-lift factor of the ogee planform is the 
lowest of the three. The highest drag factor is 
obtained with the blunt ogive planform, 20 percent 
higher than the ogee; however, since only moderate 
camber is required, that value may be more easily 
attained in practice. 

The corresponding set of pressure distribu- 
tions for flat wings of the same planform series is 
presented in figure 5. In this figure, the surface 
of the flat plate is shown in the same coordinate 
system as used in the previous figure to facilitate 
comparisons. All of the pressure distributions 
exhibit the extreme leading-edge pressures associ- 
ated with linear theory analyses of subsonic 
leading-edge wings. In the supersonic leading-edge 
portion of the blunt ogive wing, and the areas near 
the locally supersonic leading edge of the ogee 
planform, a general reduction of the pressure- 
distribution peaks may be observed. The drag-due- 
to-lift factors of the flat -plate solution repre- 
sent only the drag component caused by the rearward 
inclination of the normal-force vector at angles of 
attack, and do not include any of the theoretical 
"leading-edge- suction" force. 

The two solutions of figures 4 and 5 may be 
combined to produce a conventional force coeffi- 
cient presentation, of drag-due-to-lift and 
pitching moment versus lift, if a design lift coef- 
ficient is selected so that a particular camber 
surface may be defined. This is done in figure 6 
for three wings of the series. In this comparison, 
the wing shape in each case is that required by the 
optimum combination of loadings for minimum drag at 
a design lift coefficient of 0.10. The variation 
in drag coefficient and pitching -moment coefficient 
with increments in lift coefficient away from 0.10 
is defined by the flat -wing characteristics, using 
a subsequently discussed superposition technique. 

The lowest drag-due-to-lift factor, highest degree 
of static stability, and largest value of zero-lift 
moment are obtained with the ogee planform. (For 
this comparison, all moment coefficients were based 
on the wing length and taken about the 50-percent 
length.) However, the camber surface for the ogee 
planform also requires the most severe departures 
from the z = 0 plane, with corresponding depar- 
tures from linear theory assumptions. 

The selection of the planform series of fig- 
ures k, 5, and 6 was arbitrary; the use of other 
selection criteria would likely result in different 
conclusions regarding the planforms. The illustra- 
tion intended is that the numerical methods offer a 
useful tool in studying tradeoffs between various 
planform geometries established by a set of basic 
ground rules. 

The use of the numerical solutions is not 
restricted to the design application, where an 


optimum or other specified loading is desired. The 
general case of the inverse solution may also be 
handled, that of obtaining the pressure distribu- 
tions on a specified warped wing surface throughout 
a lift-coefficient range. Such a capability is of 
great utility, both with respect to obtaining the 
theoretical aerodynamic characteristics of a design 
surface at off -design Mach numbers, and with 
respect to estimating the effect of modifying a 
theoretical surface geometry because of imposed 
physical constraints. 

The construction of a drag polar for a warped 
surface utilizes a superposition technique, 
involving the calculation of the drag and lift 
coefficients of the surface at a selected angle of 
attack; the calculation of the drag and lift coef- 
ficients of the corresponding flat-wing planform 
per unit angle of attack; the interference drag of 
the flat-wing pressures acting on the warped sur- 
face; and the interference drag of the warped-wing 
pressures acting on the flat-wing surface, all of 
which are readily handled by the numerical approach 
The technique of combining these force coefficients 
into a composite drag polar for the warped surface 
is illustrated in figure 7. This composite polar 
considers drag-due-to-lift only; however, with the 
addition of appropriate drag increments for thick- 
ness effects and skin friction, a complete polar 
can be constructed for a wing or wing body. 

The construction of the Cm - curve for 
the warped wing is more direct, consisting of A 
shifting the flat-wing curve to pass 

through the moment coefficient of the warped wing 
at its corresponding lift coefficient. 

A comparison between numerical solution force 
coefficients and corresponding experimental data is 
presented in figure 8 for a wing-body with an | 
uncambered ogee wing. Experimental data for this 
configuration was obtained at three Mach numbers;® 
the numerical method data were calculated for the 
complete wing-body planform shown in the ( inset 
sketch. Agreement between the experimental and 
numerical method data is quite good, both in the 
moment and drag-due-to-lift comparisons. 

Similar comparisons between the numerical 
solution and experiment for both flat and warped 
wings, at and off-design Mach numbers, are pre- 
sented in figures 9 and 10 for arrow and ogee plan- 
forms. The warped wings of these figures were 
basically designed according to a restricted twist 
and camber approach,l>5 which limits the allowable 
degree of surface warpage in view of real flow con- 
siderations. However, even with these restric- 
tions, a singularity^ occurs in the camber surface 
near the wing root which must be adjusted empiri- 
cally to provide a reasonable root incidence angle. 
In calculating the theoretical aerodynamic charac- 
teristics of these warped wings the actual camber 
surface shape that was tested was used. 

The force coefficient comparison for the arrow 
wings is shown in figure 9* The wings are the flat 
and 0.08 design-lift-coefficient versions of a 
70° sweepback family. 5^6 Data are presented for 
the design Mach number 2.0 and off -design Mach num- 
bers of 1.6 and 2.2. (in this and subsequent fig- 
ures, the ACp of the experimental data is with 
respect to the minimum drag coefficient of the 
flat -wing configuration, and the ACm is with 



respect to the zero- lift pitching -moment coeffi- 
cient of the flat-wing configuration.) The numer- 
ical method predicts reasonably well the experi- 
mental drag-due-to-lift variation of the two wings 
with Mach number, particularly with respect to the 
reduction of the drag-due-to-lift difference 
between the flat and warped wings with increasing 
Mach number. Reasonably good agreement between 
experiment and the numerical solution is also 
obtained in the zero-lift moment increment AC^ 

between the two wings; however, nonlinearities in 
the experimental pitching-moment data make agreement 
poor at lift coefficients above approximately 0.10. 
It may be observed that the flat-wing drag measure- 
ments at lift coefficients on the order of 0.10 are 
lower than the theoretical estimate (Mach numbers 
1.6 and 2.0). This behavior is believed caused by 
a leading-edge vortex, 7 not included in the linear 
theory analysis. 

Also in figure 9, the drag-due-to-lift enve- 
lope polar of the three -component optimum loading 
is shown at Mach number 2.0 for the arrow wing 
planform. The increment in drag coefficient 
between the optimum envelope and the warped-wing 
polar at the design lift coefficient of 0.08 is a 
measure of the theoretical drag penalty of the 
restricted twist and camber wing, as opposed to a 
theoretically optimized shape. However, since 
experimental attempts to verify the applicability 
of the unrestricted optimum concept of linear 
theory have been largely unsuccessful,^- this theo- 
retical drag penalty is rather meaningless. Con- 
siderably more success has been obtained in exper- 
imentally verifying limited forms of the optimum 
theory; the restricted twist and camber approachl 
is one reasonably successful technique of limiting 
the local surface slopes of a wing camber surface 
so that linear theory can be used to predict the 
aerodynamic characteristics of the wing. 

The ogee planform comparison is presented in 
figure 10 with force coefficients from the numer- 
ical method and unpublished experimental data shown 
for a wing-body configuration. (The wing planform 
for the numerical solution was again taken as the 
complete projected wing -body shape.) Both flat and 
warped wings were tested, the warped version being 
designed initially for Cl = 0.10 at a Mach num- 
ber of 2.2, using the numerical method 2 and the 
restricted twist and camber approach. However, 
empirical modifications to the warped wing were 
incorporated in the course of the testing, so that 
it was necessary to scale the wing shape of the 
final test configuration from the model for the 
inverse method calculations. The numerical method 
drag polar is seen to agree reasonably well with 
the experimental data, particularly at the higher 
Mach numbers; an apparent vortex effect is evident 
in the drag data of both wings at Mach number 1.8. 

In the pitching -moment comparison, the numerical 
method is seen to predict quite well the increment 
in moment between the flat and warped wings, and to 
estimate the stability levels reasonably well. 
However, nonlinear pitching -moment effects may be 


observed in the experimental data at the higher 
lift coefficients, particularly in the case of the 
warped wing. 

Concluding Remarks 

While any technique employing linearized 
theory must be employed Judiciously in order for 
the solution to be applicable in the real flow 
case, reasonable agreement obtained between data 
from experiment and the numerical methods for var- 
ious wing configurations has verified the general 
usefulness of the numerical techniques. The use of 
the numerical methods to extend the application of 
linear theory analyses to wings having a wide vari- 
ety of planform and surface shapes offers signifi- 
cantly expanded possibilities for the development 
of aerodynamically optimized supersonic aircraft. 
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Figure 1.- Wing coordinate system. 
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Figure 2.- Numerical solution of linearized theory. 
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Figure 3.- Comparison of numerical method with superposition solution. 





1 1 — _l — I — 1 — ■ — 

0 .2 .4 .6 .8 


y 

b/2 


1.0 


0 .2 .4 .6 .8 1.0 

y 

b/2 


0 .2 .4 .6 

y 

b/2 


Figure 4.- Camber surface for a specified loading. 


.8 1.0 


FLAT PLATE 







UNCAMBERED WING-BODY 



Figure 10.- Comparison of numerical method with experiment. 



